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It is shown that the Goldstone modes associated with a broken continuous symmetry lead to 
anomalously large fluctuations of the zero field order parameter at any temperature below T c . In 
dimensions 2 < d < 4, the variance of the extensive spontaneous magnetization scales as L 4 with 
the system size L, independent of the order parameter dynamics. The anomalous scaling is a 
consequence of the l/q 4 ' d diver gence of the longitudinal susceptibility. For ground states in two 
dimensions with Goldstone modes vanishing linearly with momentum, the dynamical susceptibility 
contains a singular contribution (q 2 — c 1 ; 2 /c 2 ) _1//2 . The dynamic structure factor thus exhibits a 
critical continuum above the undamped spin wave pole, which may be detected by neutron scattering 
in the Neel-phase of 2D quantum antiferromagnets. 



It is one of the basic properties of any thermodynamic 
system that the fluctuations Var A = (A 2 ) — (A) 2 ~ V 
of an extensive variable A scale linearly with the system 
volume V. This property guarantees that the intensive 
variables A/V are self-averaging, with rms fluctuations 
vanishing proportional to V^ 1 ^ 2 . Physically, the exten- 
sive nature of Var A~ V is due to the existence of a finite, 
microscopic correlation length £. A large system can thus 
be partitioned into an extensive number of V/^ d subvol- 
umes which are statistically independent. Since the vari- 
ance in each sub volume is finite, the central limit theo- 
rem then quite generally implies a Gaussian distribution 
for A with a variance of order V, in agreement with the 
standard Einstein theory of fluctuations in macroscopic 
thermodynamics. 

The above argument indicates that the linear scal- 
ing Var A ~ V may break down only at a critical point 
of a continuous phase transition where the correlation 
length £ diverges. In fact, from the standard relation 
Var M — V ■ T x between the fluctuations of the total 
'magnetization' M and the corresponding linear suscep- 
tibility, the finite size scaling x(T c , L) ~ L 2 ~ v of the sus- 
ceptibilty right at T c implies a nontrivial dependence 
VarM ~ on S y S tena size L at the critical point. 

Below, it will be shown that anomalous fluctuations of 
the order parameter are present not only at T c but in 
fact at any temperature below T c provided the broken 
symmetry is continuous. This is a result of the pres- 
ence of Goldstone modes, which imply that correlations 
decay algebraically below T c with exponents that are in- 
dependent of temperature. A specific example of recent 
interest are fluctuations of the condensate number Nq 
in a Bose Einstein condensate. Using a weak coupling 
Bogoliubov approach, Giorgini et.al have shown Q that 
VarTVo ~ T 2 L 4 scale anomalously at low temperatures. 
In fact, this result also applies to strongly interacting 
superfluids, being essentially a consequence of the long 
wavelength phase fluctuations Q . In the present Letter, 
we will show, that the anomalous fluctuations in a Bosc- 
Einstein condensate are just one particular example of a 
rather general phenomenon which appears in any phase 



with a broken continuous symmetry. Consider, for ex- 
ample, an isotropic ferromagnet in zero field below T c . 
Because of the invariance under spin rotations, it costs 
no energy to rotate the direction of the magnetization 
vector. Provided that the microscopic interactions are 
short ranged, this implies that the transverse suscepti- 
bility Xx(?) diverges precisely as 1/q 2 at small wavevec- 
tors q —* 0. As a result, there are strong fluctuations 
of the direction of the magnetization, leading to a com- 
plete destruction of long range order at finite temperature 
in dimensions d < 2, the well known Mermin-Wagner- 
Hohenberg theorem. Regarding the magnitude of the 
magnetization, the naive expectation is, that its fluctu- 
ations are just like that of a standard thermodynamic 
variable because there is a finite restoring force for devia- 
tions from the equilibrium value. However, as was noted 
a long time ago by Patashinski and Pokrovski the 
inevitable coupling between longitudinal and transverse 
order parameter fluctuations entails that the longitudi- 
nal susceptibility is also singular at small wave vectors, 
diverging as X||(<z) ~ l/? 4 ~ d in dimensions 2 < d < 4. 
Our aim in the following is to show that 

a) the l/q 4 ~ d divergence of the longitudinal suscepti- 
bility leads to anomalous fluctuations Var M s ~ L 4 of the 
zero field order parameter in 2 < d < 4 and T ^ for an 
arbitrary broken continuous symmetry phase. For super- 
fluids, the corresponding relative condensate fluctuations 
are universal at low T. 

b) there is an analog of the singular nature of x\\ f° r zero 
temperature phases with a broken continuous symmetry. 
In two dimensions and with Goldstone modes whose fre- 
quencies ui = cq vanish linearly with momentum q = \q\, 
the dynamical susceptibility has a singular contribution 
(q 2 — cj 2 /c 2 ) -1 / 2 as noted first by Sachdev on the basis 
of a l/N expansion 0. The dynamic structure factor 
of 2D quantum antiferromagnets thus exhibits a critical 
continuum above the standard <5-function spin wave peak. 
This effect may be observed in high resolution neutron 
scattering experiments. 

As an effective description of an arbitrary phase with 



2 



a broken continuous symmetry at finite temperature we 
use the nonlinear a- model It describes directional 
fluctuations of an order parameter ip(x) — mi° Sl(x) with 
a fixed magnitude mi in terms of an TV-component unit 
vector Ct(x). At zero external field, the effective action 
for the fluctuations of ft is 



S[Sl] = ^ 

L j 2T 



d d x[Vtt{x) 



\n{xt 



1, 



(i) 



with the spin stiffness (or helicity modulus) p s as the 
single phenomenological parameter. In a finite system 
of volume V = L d , there is of course no spontaneous 
magnetization at zero field. The breaking of a continuous 
symmetry below T c however shows up in the integrated 
0(-/V)-symmetric correlation function at finite volume 



d d x $(x) ■ $0)) L = Vm 2 L -> Vm 



(2) 



defining an intensive nonzero order parameter m 2 L which 
approaches the spontaneous magnetization m 2 of the 
infinite system in the thermodynamic limit L — > oo. 
For superfluids this is just the number of particles in 
the condensate. In dimensions d > 2, where m s is 
nonzero at finite T, the leading long distance behaviour 
of the two point function G(r) = (ip(x) ■ ip(0)) may 
be obtained from a simple Gaussian spin wave calcula- 
tion. Following a standard procedure |y] , the unit vector 



Q(x) — I n(af), Y 1 — Tl(x) 2 J is decomposed into N — I 
'transverse' Goldstone fields H a (x) , a = 1, . . . N — 1 and 



a longitudinal component Q,n(x) = y 1 — II(:r) 2 . At low 
enough temperatures, where spin wave interactions can 
be neglected, the II fields are Gaussian random variables 
with variance (IL a (q)Il b (q')) = 8 a ^5 q - q >T/p s q 2 . The 
zero field correlation function below T c 



G(r) =ml[l + C\\ (r) + (N - I) • Cj_(r) 



(3) 



is thus split into longitudinal and transverse parts, with 
m 2 = G(oo) the renormalized value of the spontaneous 
magnetization. To lowest nontrivial order in the small 
fluctuations II, the transverse correlation function decays 
proportional to T/p s r d ~ 2 , as expected from the standard 
1 /q 2 divergence of the transverse susceptibility 



X ±(q)=m 2 s C ± (q)/T 



p s q 



(4) 



in the symmetry broken phase. In leading order, the lon- 
gitudinal function is given by (c stands for 'connected') 



q,( r )«i(iP(x)fl a (o)) c = ^ 



C 2 ± (r) (5) 



and thus again decays algebraically as \/r 2 ( d ~ 2 \ Con- 
trary to the naive mean field picture, where the londitu- 
dinal susceptibility X\\(l) below T c is finite as q — > 0, the 



slow decay of C\\(r) implies that X\\(q ~ * 0) ~ T j p 2 s q A ~ d 
is divergent in d < 4 4]. While the relation (5) is valid 
only to leading order in an expansion in powers of II, the 
behavior Cu (r) ~ r -2(d-2) - 1S - m ^ &c ^. ex p ec ted to be exact 
at arbitrary temperatures below T c [(| , consistent with a 
rigorous correlation inequality for N = 2 0- 

To discuss the fluctuations of the spontaneous mag- 
netization in a finite system at zero field, we define a 
fluctuating, extensive variable (dl — d d x\) 



M, 



1 

V 



dl / d2 ip(l) ■ V>(2) 



(6) 



with thermal average (M s 
then determined by the connected four point function 



Vm L . Its fluctuations are 



U4(1234) = (V»(l)-V(2)V(3)-^(4))c 



(7) 



Expanding this consistently up to order II 4 , one finds 
that in an infinite system this function may be expressed 
simply as 

U4 (1234) = m 4 ^— 

■ [C x (ri 3 ) - C ± (r u ) - C ± (r2 3 ) + C ± (r 24 )} 2 . (8) 

In order to obtain the scaling of Var M s in a finite system, 
we switch to a momentum representation and replace in- 
tegrals J q by discrete sums V~ x ' over wavevectors. 
The g = contribution is excluded since it describes an 
irrelevant global rotation of the order parameter ip(x). It 
is then straightforward to show that 



VarM s = 2m 4 s (N -I) 



(9) 



which is proportional to L 4 in 2 < d < 4 by simple di- 
mensional analysis j3| • Defining a numerical coefficient B 
by J2 q V 4 = 5L 4 /8tt 2 , we find that B = 8E 3 {2)/tt 2 = 
0.501 for a 3D box with Dirichlet boundary conditions. 
Here E d {t) = E^=i,... n<J =i (»? + ... + n 2 )- 1 is the gen- 
eralized Epstein zeta function, convergent for d < 2t (for 
periodic boundary conditions B = 0.8375 |3j). 

Although our derivation of the general result (9) is 
based on an expansion in powers of II and thus appears 
to be restricted to low temperature, the exponent in Var 
M s r~j L 4 is universal below T c just like the gr _ ( 4_£i ) di- 
vergence of the longitudinal susceptibility. Similarly, the 
temperature dependence of Var M s will be (T/p s (T)) 2 for 
arbitrary temperatures below T c , vanishing proportional 
to T 2 at low temperatures since p s (T = 0) is finite. This 
follows from the fact that at any temperature T in the 
broken symmetry phase, the dominant finite size depen- 
dence is determined by the leading low energy constant 
in the effective field theory for fluctuations of the order 
parameter, which is precisely p s as defined in equation 
(1). New effective constants only enter in higher order 



3 



terms, as discussed for instance in the context of chiral 
perturbation theory in QCD From eqn.(9), it fol- 

lows that the relative fluctuations Var M S /(M S ) 2 scale 
as £~ 2 ( d_2 ). The spontaneous magnetization is there- 
fore self-averaging in d > 2 although weaker than ex- 
pected naively, unless d > 4. In this context, it should be 
mentioned that, in the case of a non-Abelian symmetry 
N > 3, the spin stiffness p s has anomalous fluctuations 



(Ps) 2 



(N 



T 



(10) 



on its own, as shown by Chakravarty |ll| . They are very 
similar to the fluctuations of the spontaneous magnetiza- 
tion discussed here and indeed arise from the same kind 
of Goldstone anomalies. A special situation appears in 
homogeneous superfluids (N = 2) at low temperatures. 
Since translation invariance requires the superfluid den- 
sity to be equal to the full density n, the associated stiff- 
ness p s (T — > 0) = ft n/m is independent of the interac- 
tion. Equation (9) then leads to the remarkable result 
that the relative fluctuations of the condensate number 
at low temperature (d = 3) 



,. Var No 
hm — - 

T-0 (AT )2 



D 



(nX 2 T LY 



(11) 



are completely universal, depending only on density, sys- 
tem size and the thermal wavelength At — h/y/2irmT, 
but not on whether the superfluid is weakly or strongly 
interacting Finally, it is important to realize that 

- as in the case of the Mermin-Wagner-Hohenberg theo- 
rem - the result (9) is independent of the order parameter 
dynamics. It applies equally to say ferro- or antiferro- 
magnets even though the temperature dependence of the 
average order parameter (M s ) is, of course, very different 
in both cases. This may be shown either by a microscopic 
derivation of equation (9) on the basis of noninteracting 
spin wave theory around a perfectly ordered ground state 
or, alternatively, by using a quantum mechanical gen- 
eralization of the nonlinear a— model which describes 
magnons with dispersion either u> ~ q 2 or u> ~ q. 

As shown above, the anomalous fluctuations (9) of the 
order parameter are a consequence of the q~( 4 ~ rf ) diver- 
gence of the longitudinal susceptibility in d < 4. In 
the following we want to discuss the analog of this phe- 
nomenon in zero temperature phases where a continuous 
symmetry is broken. In this case the specific dynamics 
of the order parameter is important. Here, we assume 
that the Goldstone modes have a linear spectrum u> = cq 
as in superfluids and antiferromagnets. An effective de- 
scription of the ordered phase can then be obtained from 
a quantum mechanical nonlinear er— model, with a unit 
vector Q,(x, r) which depends both on space x and imagi- 



nary time r e [0, j3K\. The corresponding effective action 



sin] 



2ft 



an 



dr 



d d :i 



(12) 

has the spin wave velocity c as the only additional pa- 
rameter. Together with the renormalized value m 2 a of 
the long range order, p s and c completely determine 
the low energy properties of the ordered phase. As an 
example, this model applies both to the Neel ordered 
state of 2D quantum antiferromagnets discussed exten- 
sively in the context of high temperature superconductors 
[l3| or to the superfluid phase of cold atoms in an opti- 
cal lattice which has been realized recently in 3D |l4| . 
Using again a lowest order expansion in powers of the 
small fluctuations II(x, r) in the standard decomposition 

n = (n, a/i 

at T — in two dimensions is given by [If 



1 — ft 2 ) , the transverse correlation function 



C±(x,t) 



he 



4np s y/r 2 + (cr) 5 



(13) 



Analytic continuation to real time t = —it and Fourier- 
transformation give the standard form of the transverse 
dynamical susceptibility in any dimension 

(14) 



Ps 



q 2 — uj 2 /c 2 

It leads to the expected quasiparticle pole at ui = cq, 
reflecting the presence of undamped antiferromagnetic 
spin waves. Considering the longitudinal correlations, 
the analog of the relation (5) again applies to leading 
order. As a result, the longitudinal dynamical suscepti- 
bility in 2D turns out to be 

(N - l)m 2 hc 1 



16p 2 



\/ q 2 — uj 2 /c 2 



(15) 



which has a branch cut rather than a simple pole. For- 
mally the result (15) is completely analogous to the 1/q 
divergence of the classical static susceptibility x\\{q) m 
d = 3 discussed above. Indeed the effective dimensional- 
ity of the 2d quantum antiferromagnet is d + z = 3 and 
the dependence on frequency w is dictated by the formal 
Lorentz-invariance of the model (12). 

In order to relate these results to directly observable 
quantities, we consider the rotationally averaged stag- 
gered susceptibility (for antiferromagnets N = 3) 

N — 1 1 
Xs(q,u) = N X±(q,u) + —X\\(<l,u). (16) 

For frequencies Tllo 3> T, its imaginary part is equal to 
the dynamic structure factor S(q,u>) at w > up to a 
factor 2ft, giving 



S{q,u)=2m 2 s t 



N — 1 

N 



-5{u-cq) 



gj 0(w- cq) 
16 



c 2 q 2 



(17) 
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The longitudinal fluctuations of the Neel order thus lead 
to a critical continuum above the spin wave pole at 
u> = cq, which decays only algebraically. The contin- 
uum results from the decay of a normally massive ampli- 
tude mode with momentum p into a pair of spin waves 
with momenta q and p — q, which is possible for any 
lu > cq, with a singular cross section because of the 
large phase space. The amplitude mode is thus com- 
pletely overdamped in two dimensions ||. The relative 
weight of these fluctuations compared to the dominant 
transverse contribution is determined by the dimension- 
less parameter with £j = hc/p s the Josephson cor- 
relation length. This length controls the crossover from 
the Goldstone regime q^j <C 1, where the spin dynam- 
ics is well described by small fluctuations around a Neel 
ordered ground state to quantum critical fluctuations at 
q£j 1. The quantum critical regime has a dynamical 
susceptibility of the form 0] 



Ps V 27r 



A 



Q 



{q 2 -^/c 2 ) 1 '^ 2 



(18) 



with an amplitude Aq close to one and an exponent 
77 which is nearly zero. Now, in the regime where the 
ground state exhibits strong Neel order, the Josephson 
correlation length is only several lattice constants 0] 
and thus the quantum critical regime is hardly accessi- 
ble. In turn, there is a rather wide Goldstone regime 
characterized by q£j <C 1 and q£(T) ^> 1, where £(T) is 
the finite 2d correlation length over which Neel order is 
lost at finite temperature. Due to the exponential depen- 
dence Hal 



£(T)~£jexp(27rp s /T) 



(19) 



the relevant range of wavevectors £ -1 (T) < q < £J is 
rather wide at low temperatures. Experimentally, spin 
waves in 2D quantum antiferromagnets have been ob- 
served by inelastic neutron scattering in the vicinity of 
the Neel ordering wave vector (tt,it). In a constant u> 
scan, the first term in (17) gives rise to a sharp peak 
at q = uj/c with an amplitude ~ 1/q ^?|]- The second 
contribution in (17) due to the longitudinal fluctuations 
implies an additional algebraic tail towards smaller wave 
vectors, provided q is in the Goldstone regime. At room 
temperature, where £(T) is several hundred A and with 
typical values of £j this gives a range 0.005 < ^(A- 1 ) < 
0.05. Given the resolution in ref. [17j, detection of this 
algebraic tail appears very difficult, however high reso- 
lution measurements may be able to observe the addi- 
tional contribution from the longitudinal spin fluctua- 
tions, which apparently behave like the critical fluctua- 
tions (18), however with a rather large exponent rj = 1. 

In summary, it has been shown that Goldstone modes 
associated with a broken continuous symmetry lead to 
fluctuations of the zero field order parameter which scale 
proportional to L A at any temperature below T c . For 2D 



ordered ground states like the Neel phase of high temper- 
ature superconductors, the underlying longitudinal fluc- 
tuations lead to an additional critical continuum beyond 
the simple spin wave pole in the dynamic structure factor, 
whose detection requires high resolution neutron scatter- 
ing experiments. 
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